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1 INTRODUCTION 



ABSTRACT 

We investigate the effects of weakly and strongly stratified turbulent convection on the mean 
effective Lorentz force, and especially on the mean effective magnetic pressure. Earlier studies 
with isotropically forced non-stratified and stratified turbulence have shown that the contribu- 
tion of the turbulence to the mean magnetic pressure is negative for mean horizontal magnetic 
fields that are smaller than the equipartition strength, so that the effective mean magnetic 
pressure that takes into account the turbulence effects, can be negative. Compared with earlier 
cases of forced turbulence with an isothermal equation of state, we find that the turbulence ef- 
fect is similar to or even stronger in the present case of turbulent convection. This is argued to 
be due to the anisotropy of turbulence in the vertical direction. Another important difference 
compared with earlier studies is the presence of an evolution equation for the specific entropy. 
Mean-field modelling with entropy evolution indicates that the negative effective magnetic 
pressure can still lead to a large-scale instability which forms local flux concentrations, even 
though the specific entropy evolution tends to have a stabilizing effect when applied to a stably 
stratified (e.g., isothermal) layer. It is argued that this large-scale instability could be important 
for the formation of solar large-scale magnetic structures such as active regions. 
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Such strong fields are also hard to produce by a turbulent dynamo 
dGuerrero & KapylifeOllI) . 



The spatial and temporal coherence of the large-scale magnetic 
field of the Sun is manifested by sunspots appearing within a cer- 
tain range of latitudes from one cycle to the next. A hydromag- 
netic dynamo is commonly held responsible for the generation 



Parker 



and maintenance of large-scale magnetic fields ( cf. Moffatt 



1978 



2004 



■ 1979t iKrause & Radlerll 198(1 : fRiidiger & Hollerbach 
Brandenburg & Subram anian 2005k Some mode l s Jparker 119551 
1982L Il984l; ISpiegel & Weissll 1980l; ISpruitl 1 198 ll : ISchussler etafl 
19941 : iDikpati & Charbonneaulll999h . known as flux transport dy- 



namos, rely on the existence of strong magnetic flux tubes at the 
ba se of the convec tion zone or somewh a t belo w, see also reviews 
by iHughej d2007t) and iTobias & Weissl feOOTt) . These concentra- 
tions of magnetic field are thought to become unstable once the 
field strength exceeds a critical value. The subsequent rise of mag- 
netic flux tubes to the surface is used to explain active regions and 
sunspots. Such models, however, face a number of serious issues: 
firstly, the require d strength of the magnetic flu x tubes is of the or- 
der of 10 5 gauss ( ID'Silva & Choudhurilll993l) . which is expected 
to be a hundred times the equipartition strength which is at odds 
with estimates that the tachocline becomes unstable already when 
fields of the order of 10' 



chocline becomes unstable already when 
gauss are present (e.g. lArlt et al1l2005h . 



An alternative scenario for the large-scale solar magnetic field 
is that it is maintained within the convection zone by a distributed 
dynamo, which generates diffuse sub-equipartition strength mag - 
netic fields (e.g. lStixll976j : lBrandenbur3l2005l:lKapvla et al.l2006h . 
Unlike the flux-transport dynamo, a distributed dynamo does not 
directly explain the existence of sunspots and active regions. The 
alternative idea that sunspots have their origin within the convec- 
tion zone and that they might thus be shallow phenomena is sup- 
ported by observations showing that the rotation rate of the Sun, 
as measured by sunspots, depends monotonically on their age so 
that young spots are the fastest and oldest spots are the slowest 
( Pulk kinen & Tuominedl9l^ :lBrandenburg l2005l) . If one imagines 
sunspots floating in the plasma, the rotation rate of the youngest 
spots corresponds to roughly that at r = O.95i?0. The decreas- 
ing rotation rate as a function of age is consistent with older spots 
being anchored at increasingly higher layers where Q is smaller 
du e to its negative radi al gradient near the surface; see Fig. 4 
of iBenevolenskava et alj Jl999h . This suggests that sunspots may 
form near the surface of the Sun rather than through the buoyant 
rise of coherent flux tubes from the tachocline. This is therefore 
compatible with the distributed dynamo picture provided the dif- 
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fuse fields within the conve ction zone can for m concentrations like 
sunspots near the surface (Brandenburg| |2005l) . 

A promising mechanism that can form strong concentra- 
tions from an i nitially unifo rm magnetic field was suggested by 
iKleeorin et all d 19891. Il990h who considered the effects of tur- 
bulence or turbulent convection on the large-scale Lorentz-force. 
This work has been elaborated upon in a number of subse- 
quent papers dKleeorin & Ro gachevskii 1994; Kl eeorin et al J 19961 ; 
iRogachevskii & Kleeorin! I2007T ). They find that, for a given range 
of large-scale magnetic field strengths, there is a negative turbu- 
lence contribution to the mean magnetic pressure, and the effective 
mean magnetic pressure that accounts for the turbulence effects can 
be negative. This results in an excitation of a large-scale instability. 
The growth rate of the instability inc reases as a function of den- 
sity stratification ( Ke rnel et alj2012bl> . In the Sun the density drops 
steeply in the outermost layers, which favours the development of 

this instability there. 

The strongly stratified large -eddy simulations of lUstvugovl 

( 2009) and lKitiashvili et alj (2010) may already have detected mag- 
netic flux concentrations in turbulent convection formed from ini- 
tially uniform vertical magnetic fields. Note also that a segrega- 
tion into strongly and weakly magnetized regions in magneto- 
convection has b een observed in numerical simulations at large 
aspect ratios by dTao et all 1 19981 ; iTian & Petrovavl l2012h . which 
may have its origin in some me an-field effect of the type con- 
sidered here. Also simulations of IStein et alj JIoT 1) with a hori- 
zontal uniform field at the bottom of the domain show emergence 
of ma gnetic flux structures, while a number of numerical studies 
(e.g.lSchussler & V6glerll2006t iMartmez et aLll2008t iRempel et al] 
|2009|) use strongly nonuniform fields as initial or boundary con- 
ditions. The origin of such nonuniform fields is therefore not ad- 
dressed in these latter studies. 

Direct nume rical simulations (DNS) of homogeneous 
2010, hereafter BKR) and density stratified 
20121 hereafter BKKR) forced turbulence have 



Brandenburg et al 



( Brandenburg et al 



shown that the effective magnetic pressure is negative for field 
strengths below about 40% of the equipartition value, provided 
the magnetic Reynolds number exceeds unity. However, definitive 
proof of an instability associated with the negative effective mag- 
netic pressure phenomenon came only more recently with DNS of 
forced turbulence that have sufficiently many t urbulent eddies in th e 
simulation domain ^Brandenburg et al.ll201ll ; iKemel et al]|2012ij) . 
This work has onl y become possible due to earlier DNS (BKR, 
IKemel et alj|2012hf) exploring first the relevant parameter regime. 
In the present study we investigate the effect of turbulent convec- 
tion on the effective mean Lorentz force in DNS and study the for- 
mation of large-scale magnetic structures in mean-field models. 



2 EFFECTIVE MEAN LORENTZ FORCE 

In this section we state the underlying equations, highlighting the 
difference to earlier work where anisotropic contributions from 
gravity were either weak or absent. 



2.1 Governing equations 

In this study we are mainly interested in the effects of turbulent 
convection on the mean Lorentz force. To this end we consider the 
momentum equation, 



d_ 

dt 



pUi 



_d_ 

dx. 



n»j + pgi 



(i) 



where g is the acceleration due to gravity, 

Uij =pU i U j +S ij (p+ \B 2 ) -BiBj -2upS ij , 



(2) 



is the momentum stress tensor, U and B are the velocity and mag- 
netic fields, p and p are the fluid pressure and density, 8ij is the 
Kronecker tensor, v is the kinematic viscosity, and 



^(diUj + djUi 



(3) 



is the trace-free rate of strain tensor. Throughout this paper we have 
adopted units where the vacuum permeability fio is set to unity, 
although we do include it in some expressions for clarity. 

Neglecting correlations between velocity and density fluctu- 
ations for low-Mach number turbulence, the averaged momentum 
equation is 



d_ 

dt 



7>u t 



JL 

dx. 



P9i, 



(4) 



where p is the mean fluid density, U is the mean fluid velocity, 



ilij = Ily + Tlij is the mean momentum stress tensor split into 
contributions resulting entirely from the mean field (indicated by 
superscript m) and those of the fluctuating field (indicated by su- 
perscript f). The tensor ITL has the same form as Eq. ((2), but all 
quantities have now attained an overbar, i.e. 



Uij ^pUiUj+5 tJ 



{p+\B 2 



- 11,11, - 2vpSij 



(5) 



where B is the mean magnetic field and p is the mean fluid pres- 
sure. The contributions, IL;,-, which result from the fluctuations of 
velocity u = U — U and magnetic fields b — B — B, are deter- 
mined by 



II 



ij — pUiUj + 2<Jij 



+ ±5yb 2 - hbj 



(6) 



This contribution, together with the one from the mean field, IL^ , 
comprises the total mean momentum tensor. The contribution from 
the fluctuating fields is split into parts that are independent of the 
mean magnetic field (which determine the turbulent viscosity and 
background turbulent pressure) and parts that do depend on the 
mean magnetic field. 

In the present study we consider turbulent convection with an 
imposed uniform horizontal magnetic field, Bo = (Bo, 0, 0), that 
is perpendicular to the direction of gravity. This modifies the stress 

tensor from IL^ to II^ 8 , so only the difference, 



All EE II 



n f '° 



(7) 



de pends on the mean magnetic fi eld B and can be parameterized 
as {Rogachevskii & Kleeorin 20q3) 



AIT-j = q B B 2 BiBj - (±q p <5y + q s fog^ B 2 , 



(8) 



where g is the vertical unit vector directed along the gravity field, 
Bj — Bj/B is the unit vector directed along the mean magnetic 
field, q B ,<lp, and q g are functions of magnetic Reynolds and Prandtl 
numbers as well as the modulus of the normalized mean field, 



f3 = B/B ccl , where B = \B\, and B cq = (pu 2 



,1/2 



(9) 



is the equipartition fiel d strength. To derive Eq. (|8) , we 
use Eqs. (A22)-(A24) of Ho gachevskii & Kleeorinl l2007h . The 
parametrization ([8} follows also from symmetry arguments which 
allow us to construct a symmetric tensor with two preferential per- 
pendicular directions along the horizontal magnetic field B and 
vertical gravity field g. Such symmetric tensor is a linear combi- 
nation of symmetric tensors 5y, BiBj and gtgj. (In the case of 
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an oblique imposed magnetic field, there would be an additional 
contribution.) 

The effective mean Lorentz force that takes into account the 
turbulent convection effects, reads: 



pK = -V, ( ±B% - BiBj + Af4 



+B V 



(1-<?p)S 2 



+ 9 



i V z (q g B 2 ) 



(10) 



The analytic expressions for the nonlinear quenching functions, 
q D (P), q s (/3), andq B (/3) for turbule nt convection have been derived 
in Rogachevskii & Kleeorin | |2007| ). Their asymptotic formulae are 
given below. For weak mean magnetic fields, 4/3 <C R^ 1 ' 4 , the 
functions q p , q s , and g g are given by 



?p(/?) 
for R, 

*(/J) 



4 / , „ 4 
- lni? m + — 

5 V 45 



il 1 



7 , 16i§ 
— - a* + 



8a, 



+ 15). &0) 



-1/4 



< 4/3 < 1 these functions are 

5|ln(4/3)| + l + 32/3 2 



16 

25 

32 
15 

8a„ 



7 

a* 

3 



3tf 2 ' 



16 £§ 
9 H P 2 ' 



|m(4/3)| + - + 12/3 2 



3B P 2 ' 



while for strong fields, 4/3 2> 1, they are 



*G8) = 



j_ri+^iui^ 

6/3 2 V H 2 )^ 80/3 ' 



37ra» 



48/3 3 160,0' 



+ 



37ra, 



40/3 3 4tf 2 ,fl 2 



Here, 7? m = foW ums /?7 is the magnetic Reynolds number based 
on the integral scale of turbulent convection, to, and the root-mean- 
square (rms) value of the velocity, u rms , r\ is the magnetic diffusion 
due to the electrical conductivity of the fluid, and H p is the density 
scale height. The parameter a, characterizes turbulent convection 
and is determined from the budget equation for the total energy, 
yielding 

a' 1 = 1 + [^{Wf + »jt(VB) 2 />] IgF, , 

where vx is the turbulent viscosity, rj t is the turbulent magnetic dif- 
fusivity, _F» = u z s' is the vertical heat flux from the background 
turbulent convection, and s' are fluctuations of the specific entropy. 



2.2 Turbulent contributions to effective Lorentz force 

To study the effects of turbulent convection on the Reynolds and 
Maxwell stresses, and hence on the effective Lorentz force from 
the mean field, we need to determine the functions q P (/3), q s (P), 
and o g (/3) in DNS. Allowing here for the possibility of small-scale 
dynamo action, Eqs. l[6](-([8]( yield 



ApUiUj + ^SijAb 2 — Abibj 

= [q s BiBj - (iq p <Jy + q s gi gj )] B 2 . 



Here, AjMiuJ = pUiUj — TTuoT/uoJ and Abibj = bibj — boiboj, 
where subscripts indicate values in the absence of the mean mag- 
netic field. To obtain 3 independent equations for the 3 unknowns, 
we multiply Eq. i ll It : 

(i) by gtgj, 

A^f+iA^-Afef^-^igp+gg^B 2 , (12) 

(ii) by BiBj (defining — u ■ Bo and b^ — b ■ Bo): 



Apu% + § Ab 2 - Ab% = - [ \q p -q a )B" 



(13) 



(iii) and compute the trace of Eq. i ll It . 
A^ + § Ab 2 " = - (§q p + g g - g s ) B 2 
Equations d 1 2b — d L4b yield the functions g p (/3), g s (/3), and q g (/3): 



(14) 



q p (P)B = -2Apu% - Ab 2 + 2A6 2 , (15) 
g s (/3)B 2 = -A~p~u> y + A^ 2 +A6|- A&f, (16) 
q s (fi)B 2 = -A^f + A^f + A&i- Abf. (17) 
Using Eqs. 

we determine the functions q P (/3), q s (/3), 
and ffg(/3) from DNS in Sect. [3] In all those cases, overbars denote 
averages over x, y, and t during the statistically steady state, corre- 
sponding to a time interval At of up to a thousand turnover times, 
i.e., Atu tms kf = 0(1000). 



3 DIRECT NUMERICAL SIMULATIONS 
3.1 DNS model 

We use two different setups when studying the effects of convection 
on the mea n Lorentz force: (i) a weakly stratified model, similar to 
that used in iKapyla et alj 1201CI) . without ove rshoot layers, and (ii ) 
a strongly stratified setup, similar to that in IKapyla et al] ^2009), 
including upper and lower overshoot layers. In both cases we use 
a Cartesian domain with L x — L y — 5d, where d is the depth 
of the convectively unstable layer. The convective layer is situated 
between < z < d in both setups. In the strongly stratified setup 
the z coordinate runs from — 0.85d < z < 1.15. We solve a set of 
hydromagnetic equations 



OA 

— = U x B - r/poJ , 

2^ = --\?p + g + -J x B + -V ■ 2vpS, 
Dt p p p 



Dlnp 
Dt 

Dt 



= -V-U 
1 

P 



VKVT + 2vS 2 + ^J 2 



(18) 



(19) 



(20) 



(21) 



(11) 



where D/Df = d/dt + U ■ V is the advective time derivative, A 
is the magnetic vector potential, B = VxA + Bo the magnetic 
field, B = (Bo, 0,0) is the imposed field, J — Pg 1 ^ x B is 
the current density, 77 and v are the magnetic diffusivity and kine- 
matic viscosity, respectively, K is the heat conductivity, p is the 
density, U is the velocity, s is the specific entropy, and g = —gz 
is the gravitational acceleration. The fluid obeys an ideal gas law 
p = pe("i — 1), where p and e are the pressure and internal energy, 
respectively, and 7 = cp/cv = 5/3 is the ratio of specific heats 
at constant pressure and volume, respectively. The specific internal 
energy per unit mass is related to the temperature via e = cvT, 
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and the rate of strain tensor S is given by Eq. l[3}. The stratification 
in the hydrostatic initial state can be described by a polytrope with 
index m — 1 in the weakly stratified case, and stratified three-layer 
setup is described by polytropic indices (mi, 7712, mz) = (3, 1, 1). 
In the latter setup, a cooling term r coo i operates in the region 
z > 1, keeping the layer isothermal. The density contrast across 
the full domain, Ap = pt,ot/pto P , is 1.2 (runs A) and 320 (runs 
B, see Table [TJ. In the latter case, the density changes by a factor 
of roughly 13 within the convectively unstable layer. Typical flow 
patterns from both setups are shown in Fig.Q] In all simulations we 
use the PENCIL CODlfl 

The horizontal boundaries are periodic. In the weakly strat- 
ified case we keep the temperature fixed at the top and bottom 
boundaries, whereas in the more strongly stratified setup the up- 
per boundary is isothermal and a constant flux of energy is applied 
at the lower boundary by fixing the temperature gradient. For the 
velocity we apply impenetrable, stress-free conditions according to 



d z U x = d M U v = U Z = 0. 

For the magnetic field we use vertical field conditions 



B x — Bo = B v 



0. 



(22) 



(23) 



Dimensionless quantities are obtained by setting 

d = g = Po = cp = po = 1 , (24) 

where po is the fluid density at z m = ~d. The units of length, time, 
velocity, density, specific entropy, and magnetic field are then 



(25) 



d, [t] = Vd/g, [U] = Vdg 



[p] = Po , [s] = c P , [B] = \fdgp 



■opo 



The simulations are controlled by the following dimensionless pa- 
rameters: thermal and magnetic diffusion in comparison to viscos- 
ity are measured by the Prandtl numbers 



Pr=A 



Pm=^, 



(26) 
Xo V 

where xo = K/(cppo) is the reference value of the thermal dif- 
fusion coefficient measured in the middle of the layer (z m ) of the 
non-convecting hydrostatic reference initial state. The efficiency of 
convection is characterized by the Rayleigh number 

v%o \ cp dz, 

which is again determined from the initial non-convecting state at 
z m . The entropy gradient can be presented in terms of logarithmic 
temperature gradients 

V- V ad 



(27) 



1 ds\ 
cp dz J - 



(28) 

1 /7, and Hp being the 



with V = (<91nT/dlnp) Zm , V a d = 
pressure scale height at z — z m . 

The effects of viscosity and magnetic diffusion are quantified 
respectively by the fluid and magnetic Reynolds numbers 



Re 



vk 



Rm = 



Tjkf 



= Pm Re, 



(29) 



where u rms is the root-mean-square (rms) value of the velocity and 
kf — 2n/d is the wavenumber corresponding to the depth of the 
convectively unstable layer. Again, it is convenient to measure the 
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Table 1. Summary of the runs. Here, Ma 

— tirais/ Y d>g, and the imposed 
field in normalized form is given by Bo = Bo/B eq , where B cq is the 
volume averaged equipartition field. The last column gives the density con- 
trast, Ap, over the entire domain. The Prandtl number is equal to unity in 
all runs. We use grid resolutions 128 2 X 64 (Runs A1-A24), 256 2 X 128 
(Runs A25-A29), and 256 2 X 192 (Runs B1-B8). 



Run 


Ma 


Ka 


Rm 


Pm 


Bo 


Ap 


AOh 


0.068 


10 6 


— 


— 


— 


1.2 


Al 


0.068 


10 6 


11 


0.2 


0.07 


1.2 


A2 


0.066 


10 6 


11 


0.2 


0.15 


1.2 


A3 


0.059 


10 6 


9 


0.2 


0.34 


1.2 


A4 


0.055 


10° 


9 


0.2 


0.54 


1.2 


A5 


0.052 


10° 


8 


0.2 


0.78 


1.2 


A6 


0.053 


10 6 


8 


0.2 


0.95 


1.2 


A7 


0.059 


10° 


9 


0.2 


1.17 


1.2 


A8 


0.066 


10 6 


1 1 


0.2 


1.52 


1.2 


A9 


0.068 


10 6 


22 


0.4 


0.07 


1.2 


A10 


0.065 


10° 


21 


0.4 


0.15 


1.2 


All 


0.055 


10° 


18 


0.4 


0.36 


1.2 


A12 


0.051 


10 6 


16 


0.4 


0.59 


1.2 


A13 


0.048 


10 6 


15 


0.4 


0.83 


1.2 


A14 


0.059 


10 6 


16 


0.4 


1.01 


1.2 


A15 


0.061 


10° 


19 


0.4 


1.15 


1.2 


A16 


0.066 


10° 


21 


0.4 


1.51 


1.2 


A17 


0.065 


10 6 


51 


1.0 


0.08 


1.2 


A18 


0.061 


10 8 


49 


1.0 


0.16 


1.2 


A19 


051 


1Q 6 


4 1 


1.0 


39 


1 2 


A20 


047 


10 6 


37 


1.0 


0.64 


1 2 


A21 


0.047 


10 6 


37 


1.0 


0.86 


1 .2 


A22 


050 


10 6 


39 


1 


1 01 


j 2 


A23 


059 


10 6 


47 


1 


1 18 


j 2 


A24 


0.069 


1Q 6 


55 


1.0 


1.46 


1 .2 


A25h 


0.058 


4.2 ■ 10 6 


_ 


_ 


_ 


1.2 


A25 


0.058 


4.2 ■ 10 6 


92 


1.0 


0.09 


1.2 


A26 


0.044 


4.2 ■ 10 6 


70 


1.0 


0.46 


1.2 


All 


0.040 


4.2 ■ 10 6 


63 


1.0 


0.76 


1.2 


A28 


0.042 


4.2 ■ 10 6 


66 


1.0 


0.96 


1.2 


A29 


0.047 


4.2 ■ 10 6 


75 


1.0 


1.07 


1.2 


BOh 


0.032 


1.2 ■ 10 7 








320 


Bl 


0.032 


1.2 ■ 10 7 


51 


1.0 


0.01 


320 


B2 


0.031 


1.2 ■ 10 7 


49 


1.0 


0.07 


320 


B3 


0.028 


1.2 ■ 10 7 


45 


1.0 


0.23 


320 


B4 


0.027 


1.2 ■ 10 7 


43 


1.0 


0.32 


320 


B5 


0.027 


1.2 ■ 10 7 


43 


1.0 


0.41 


320 


B6 


0.026 


1.2 ■ 10 7 


42 


1.0 


0.61 


320 


B7 


0.025 


1.2 ■ 10 7 


40 


1.0 


0.78 


320 


B8 


0.025 


1.2 ■ 10 7 


40 


1.0 


0.92 


320 



magnetic field strength in terms of the equipartition value. The val- 
ues of these parameters used in different runs are given in Table [TJ 



3.2 Effective mean Lorentz force from DNS 

We now turn to DNS models of turbulent convection to determine 
the coefficients q p , q s , and q g using Eqs. dl5t-J17t. First, we per- 
form purely hydrodynamical simulations to determine the turbu- 
lent background velocity u 2 in the absence of magnetic fields. No 
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Run BOh 

i- 500.0 



Run A35h 
t= 500.0 




EV(<*C/) 1/2 
— 0.05 



Figure 1. Velocity component U z at the periphery of the domain from hydrodynamical runs A25h (Ap = 1.2, left) and BOh (Ap : 
the top and bottom slices show Uz near the top and bottom of the convectively unstable layer, respectively. 



: 320, right). In both cases 



dynamo action occurs in Runs A1-A16 and B1-B8, whereas in 
Runs A17-A24 the dynamo is growing very slowly, and in A25- 
A29 a small-scale dynamo is operating. We find that the critical 
Rm for Pm = 1 is between 50 and 60, which is almost two time s 
higher than in the case Pm = 5 studied earlier dKiipvla et al.l2008h . 
The hydrodynamical simulations have been run sufficiently long 
(roughly 250 turnover times) so that it is thermally relaxed and the 
turbulence is statistically steady. The last snapshot of the hydro- 
dynamical run is used as the initial condition for all subsequent 
simulations where a uniform magnetic field is imposed. 

Next, we consider an imposed horizontal field, Bo = 
(Bo, 0,0). Earlier numerical studies have shown that the effec- 
tive mean magnetic pressure that takes into account turbulence ef- 
fects, is negative when the mean magnetic fields are smaller than 
the equipartition strength in non-stratified (BKR) and stratified 
(BKKR) forced turbulence. In the present study we investigate this 
issue for turbulent convection with weak and strong density strat- 
ification. Let us define the dimensionless effective mean magnetic 
pressure V e s as 



Pcff = i(l-g P )B/B 2 



(30) 



where q p = q D (B / B ea ) and B e a = B eq (z). Following earlier 
work (BKKR. lKemel et alll2012a ). we characterize our numerical 
results for q p by a fit of the form 



<?p = 



q P o 



1 + B^/Bl : 



(3D 



which are determined 



where q p o and B p are fit parameters, 
by matching the shape of V c s near its minimum and where 

2 

dVcff/dB < 0, which is the range relevant to the negative mag- 
netic pressure instability (Rogachevski i & Kleeorinf2007l BKKR). 
We use the same ansatz also for q s and q s , and define in this way 
the fit parameters q s o, <?g0, B s , and B g . Note that, due to turbu- 
lent pumping effects, the mean magnetic field B in general also 
de pends on height — ev en in the absence of a large-scale dynamo; 
see lKapyla etall 1201 fj) and BKKR. 



j* 0.0 




-0.2 



B/B K 



Figure 2. Effective magnetic pressure as a function of the mean magnetic 
field from weakly stratified Runs A1-A29 with an imposed horizontal field 
Bo = Boa;. The black stars, red diamonds, blue crosses, and yellow tri- 
angles denote simulations with Rm 10, 20, 50, and 70, respectively. We 
omit points near the boundaries at z/d < 0.35 and z/d > 0.65. The 
two curves correspond to approximate fits determined by Eq. ( I3H . with 
9p0 = 40 and B p = 0.1B cq (upper curve, small Rm), and q p o = 130 and 
B p = 0.08-Bcq (lower curve, larger Rm), respectively. 



Results for the effective mean magnetic pressure Vcs from 
weakly and strongly stratified runs are shown in Figs. [2}|3] Since 
B and B cq are functions of z, we obtain for each combination of 
imposed field strength and Rm a family of solutions for q p , q s , and 
q g . We neglect points near the top and bottom boundaries to avoid 
boundary effects. The result is shown in Fig. [2] For weak stratifica- 
tion (Runs A1-A29), a negative contribution of turbulent convec- 
tion to the mean magnetic pressure is found if the mean magnetic 
field is smaller than the equipartition value. The maximum of this 
contribution is attained near B 
what stronger for larger magnetic Reynolds number. 
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Figure 3. Same as Fig.[2]but for Runs B1-B8 for Rm = 40-50. The solid 
line corresponds to a fit with q p o = 95 and B p = 0.04£? C q. 



Figure 4. q g as a function of mean magnetic field from Runs A1-A29. The 
symbols and colours are the same as in Fig. [2] The solid line applies to the 
fit parameters q g o = 200 and B s = 0.025B cq . 



In the strongly stratified runs (see Fig. |3} we also find a neg- 
ative contribution of turbulent convection to the mean magnetic 
pressure, but it is constrained to somewhat lower values (B < 
0AB cq x) of the mean magnetic field, and the effective mean mag- 
netic pressure, "P e ff, has a weaker minimum than in the weakly 
stratified case. It appears that we find universal scaling for Vets 
as a function of B/B ecl as was obtained in BKKR for stratified 
forced turbulence, provided that only data points near the middle 
(0.35 < z/d < 0.65) of the convectively unstable layer are used. 
Furthermore, our highly stratified simulations show that the mini- 
mum of V c ff and the range of the mean magnetic field in which V c s 
is negative, are roughly consistent with those found by BKKR. 

Due to anisotropy of turbulent convection there is an signifi- 
cant contribution to the effective mean magnetic pressure charac- 
terized by the term q s . This function affects the vertical component 
of the effective mean Lorentz force; see the second term on the 
right-hand-side of Eq. dlOt . Our DNS show that the dimensionless 
quantity g g is mostly positive (see Figs. [4] and [3J, which implies 
that this effect increases the negative contribution of anisotropic 
turbulent convection to the effective mean magnetic pressure. Note 
that the DNS in stratified forced turbulence of BKKR has not found 
strong anisotropic contributions as characterized by the term q g . On 
the other hand, the negative contribution of anisotropic turbulent 
convection to the effective mean magnetic tension, as characterized 
by positive values of q s , has neither been found in our DNS (see 
Figs. [6] and [7} nor in those of BKKR. 

In Figure [8] we show the magnetic field component B x from 
Run B3 with an imposed horizontal magnetic field B « 0.23B aq x. 
The structure of the magnetic field, however, does not show clear 
signs of magnetic flux concentrations in the DNS. Even after ad- 
ditional averaging over time and along the x direction no spa- 
tial modulation of the m agnetic field is seen. The simulations of 
iBrandenburg et ail (201 l]) strongly suggest that the reason for this 
is related to lack of scale separation. As demonstrated in Fig. 17 
of BKKR, at larger scale separation the turbulent diffusivity on the 
scale of the domain becomes weak enough to allow for the devel- 
opment of large-scale magnetic structures. 




1.00 



Figure 5. Same as Fig.[4]t>ut for Runs B1-B8. The fit parameters are q s o 
3andB g = 0.5/3 cq . 
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Figure 6. Effective magnetic tension parameter as a function of the mean 
magnetic field from Runs A1-A29. The symbols and colours are the same 
as in Fig. [2] The solid line applies to the fit parameters q B o = —30 and 
B s = 0.07B CQ . 
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Figure 7. Same as Fig.|6]but for Runs B1-B8. The solid line applies to the 
fit parameters q B o = —8 and B B = 0.09B cq . 



t= 2000.0 




Figure 8. Magnetic field component B x from Run B3 with imposed hor- 
izontal field Box at Rm = 45, Pm = 1, a density contrast of 320, and 
Bo/Beq = 0.23. 



4 ENTROPY EVOLUTION IN MEAN-FIELD MODELS 
4.1 Mean-field equations 

We now apply a mean-field model similar to that of BKR for adi - 
abatic stratification and those of BKKR and I Kernel et ID {2012b) 
for isothermal stratification. Both types of models are in principle 
able to display a large-scale instability provided the domain is big 
compared with the typical size of turbulent eddies, i.e., the scale 
separation ratio is large. In the mean-field calculations of BKKR it 
was shown that in models with too low scale separation ratio the 
turbulent magnetic diffusivity and turbulent viscosity (which are 
proportional to the scale of the energy-carrying turbulent eddies) 
was too large, so the instability is too weak or not excited. Even 
if the scale separation ratio is big enough, the instability can only 
develop if dV c s /d/3 2 , taken at the value of the imposed field, is 
negative inside the domain dKemel et alj|2012bl , BKKR). If these 
conditions are satisfied, the maximum growth rate of the instability 
was shown to be independent of the strength of the imposed field 
for models with isothermal background stratification. 

Whenever the instability is possible, its nonlinear develop- 
ment appears to be rather similar f or isothermal and adia batic back- 
ground stratification. In particular. iKemel et al.l 112012b) found that 
for q s — 0, the eigenmode shows no variation along the direc- 
tion of the applied magnetic field. Conversely, if the model is two- 
dimensional with no extent in the y direction, which will be assume 



here, the results are independent of the value of q a o, so we take in 
the following q B o — 0. Furthermore, for large magne tic Reynolds 
numbe r simulations at different scale separation ratios (IKemel et al.l 
l2012al) suggest q p o = 40 and (3 P = 0.05, whi ch were there f ore also 
the fiducial parameters used in the study of IKemel et al. I d2012bh 
and will therefore also be used here. The results of the DNS pre- 
sented here suggest somewhat larger values of q p o of 130 for Set A 
and 95 for Set B, but this could be a consequence of intermedi- 
ate magnetic Reynold s numbers for whic h q g o is known to reach 
a peak (see Fig. 10 of Kern el et"aill2012j ). The dependence on the 
parameter q g o has not yet previously been determined, so this will 
be done at the end of Sect. 14.31 In all other cases we keep q e o = 0. 
The imposed field strength Bo is chosen such that the minimum of 
V e f{ occurs near the top boundary and thus dP or f /d/3 2 < in the 
domain. We express Bo in units of B cq o = Bcq(0), which is the 
equipartition field strength at z = 0. 

The novel aspect of the present work is that an evolution equa- 
tion for the mean specific entropy is included. Thus, we solve the 
following system of equations for the mean vector potential A, the 
mean velocity U, the mean density p, and the mean specific entropy 
s, in the form 



— = U X B - TjTfJ-oJ, 



-UVU- =Vp + g + T M + Tf ot , 



OA 

dt 

dU 

dt 

dp 
dt 

ds 
dt ' 



-U • Vp-pV -u, 



(32) 
(33) 
(34) 



U ■ Vs - -L V ■ F + 2v T S 2 + 55^° J 2 - Ir-cooi, (35) 
pT P T 



where B — Bo + V x A is the mean magnetic field including the 
imposed field, rpr = ??t + V and iat- = v t + v are total (turbulent 
and microphysical) magnetic diffusivity and viscosity, respectively, 
the effective mean Lorentz force is given by Eq. dlOt , and the total 
viscous force, T^t ot = (2/p)V ■ (pzApS). We assume Vt/rjt = 1 
for the turbulent magnetic Prandtl number. The mean temperature 
obeys (7 — l)c p T = 7p/p = c 2 . The boundary conditions are 
stress-free for the velocities, and perfect conductor boundary con- 
ditions for the magnetic field as in previous mean field models. In 
the following we consider two types of background stratification: 
isothermal and adiabatic. 



4.2 Isothermal background stratification 

We begin by assessing the effects of e ntropy evolution in the 
isothermal models studied by BKKR and lKemel et"al] f2012tj ). In 
general, the gas will not stay isothermal, because the temperature 
changes due to adiabatic expansion and compression. Indeed, in an 
isothermally stratified layer a rising blob cools adiabatic ally, be- 
comes denser or heavier, and thus experiences a restoring force 
with the Brunt- Vaisala frequency N, where iV 2 = g ■ Vs/c p — 
(7 ~ ^)g/"/H p , with Hp being the density scale height. It turns out 
that in such a case the negative effective magnetic pressure instabil- 
ity can be stabilized. To study this in more detail, we allow for cool- 
ing term of the form r coo i = (T — To)/c p rr,whereris the mean 
temperature, To is the reference temperature of the layer, and r is a 
cooling time. For the energy flux F we assume F = —K'VT. 

In Fig. i we show the evolution of the rms value of the mean 
flow ([/ rms ) normalized by u ms — Skfr/t using kf = 2ir/H p , 
for different cooling times r and Bo = 0.1B O q0- The instability 
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Figure 9. Evolution of U Tms /u Inls for different cooling times r normal- 
ized by the Brunt-Vaisala frequency N. 



is found to operate only when rN is less than a critical value of 
order unity. The largest growth rate seen in Fig. |9]is ~ 60j?t&i- The 
earlier results with an isothermal equation of state are recovered in 
the limit tN —¥ 0, in which case a gro wth rate of « 110rj t ki is 
found; see Fig. 4 of lKemel et alj (20 12th . where the growth rate is 



normalized by (ut + rjt)k 1 . 



4.3 Adiabatic background stratification 

Owing to the stabilizing properties of stable stratification, we study 
now the evolution of the instability in an adiabatically stratified 
layer, where this stabilizing effect is absent and the squared sound 
speed is given by = g(z — Zoo). Our setup is similar to that 
of BKR, who considered a reference height z = at which ini- 
tially c s = c s o and p — p , where c s o and po are normalization 
constants. Length is normalized with respect to the density scale 
height H p o = c^/g at z = 0. This implies that = (3/2) Hpo. 
In BKR the domain extended in the x direction from —5H p o to 
+5H p o, but the resulting horizontal wavelength of the fastest grow- 
ing eigenfunction was then about half the x extent. Therefore we 
consider here a smaller domain, with —3H p o < x < 3H p o and 
—5H p o < z < H p o. As in BKR, we choose Bo/_B cq o = 0.01. 
Here we use F = — Xt'pT'Vs for the energy flux which is appro- 
priate for a turbulent layer, and v x /xt = 1 is the turbulent Prandtl 
number, where \t is the turbulent heat conductivity. 

It turns out that heating is weak, so we ignore the cooling term, 
i.e., t — y oo. In Fig.[l0]we show velocity vectors together with B y 
as well as s for three different times close to saturation. Note that 
there is a weak enhancement of s at the location where the instabil- 
ity develops a positive maximum. This mean entropy enhancement 
is associated with turbulent viscous heating, in particular the con- 
tribution ~ ^t(V • U), which is important near the surface, even 
though the magnetic flux concentration later descends to greater 
depths. 

Finally, we use this model to assess the dependence on the 
parameter g g o, which can directly contribute to the negative effec- 
tive magnetic pressure instability. According to the DNS we have 
q g0 » 200 with /3 g = Bg/B cq = 0.025 for Set A (Fig.gJ and 
g g o ~ 3 with /3 g = 0.5 for Set B (Fig.|5j. Figure [TT]indicates that 
the q s effect would be detrimental to the instability for Set A, and 
negligible for Set B. Furthermore, near q g o = the dependence of 
the growth rate on g g o is not monotonous: for g g o = 3 the growth 
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Figure 10. Velocity vectors superimposed on color scale representations 
of By (left) as well as color scale representations of s (right) for three 
different times close to saturation. Specific entropy is shown in units of 
so = 10 -4 c p . 
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Figure 11. Evolution of C/rms/ttrms for adiabatic stratification using dif- 
ferent combinations of q g o and /3g = B g /B cq . 



rate is slightly enhanced and for g g o = 10 it is decreased, but en- 
hanced for g g o = —10 by a similar amount. The saturation level is 
only weakly affected by the value of g g p. We als o checked that, as 
expected from earlier work dKemel et al. 2012b), the value of g s o 
affects neither the growth rate nor the saturation value of (7 rm s • 



5 CONCLUSIONS 

The present simulations have demonstrated that for weak strati- 
fication, and magnetic fields less than the equipartition value, a 
destabilising contribution to the mean Lorentz force is obtained 
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in the presence of turbulent convection. A similar effect is found 
for strong density stratification, although the effect is weaker and 
limited to a narrower range in magnetic fields. Our DNS results 
agree at least qualitatively with those from non-stratified (BKR) 
and strat ified (BKKR) forced turbulence and with theoretical pre- 
dictions (Rogac hevskii & Kleeorinl [20071) . although the minimum 
effective magnetic pressure can now be even more negative and the 
range where it is negative extends now to nearly 0.5-B eq . 

Such negative contributions to the effective magnetic pressure 
facilitates an instability that can lead to the generation of flux con- 
centrations from an initially uniform magnetic field. This could ex- 
plain the origin of active regions and sunspots. However, no clear 
signs of instability are found from DNS in turbulent convection. A 
possible reason is that the scale s eparation in the DNS is in suffi- 
cient, as has been demonstrated bv lBrandenburg et al.ld201ll) . who 
found conclusive evidence of the operation of the negative effec- 
tive magnetic pressure instability in DNS in forced turbulence with 
a scale separation ratio of 15, using also spatial-temporal averag- 
ing. For a scale separation ratio of 30, the effect is stro nger and flux 
concen trations can already be seen without averaging (Ke rnel et al.l 
l2012ah . On the other hand, if the scale separation ratio is as low as 
5, no flux concentrations have been found BKKR. 

So far, all studies of the negative effective magnetic pressure 
effect have only considered the case of a horizontal mean field. 
While this is the most relevant case in view of applications to stars 
with differential rotation, the case of a vertical field is also of in- 
terest and might lead to additional effects. Another important ex- 
tension of our work is to the case where small-scale dynamo action 
becomes important. This effect is expected to lower the relative 
importance of the negative effective magnetic pressure effect, al- 
though this will depend on the value of the magnetic Prandtl num- 
ber, which is yet another important parameter whose effect on the 
instability is not yet sufficiently well understood. Finally, we have 
mentioned the possibility of finite scale separation effects, which 
would mean that the scale of the magnetic structure will be impor- 
tant in determining the efficiency of the negative effective magnetic 
pressure effect. The hope is that such considerations will provide 
some insight as to why we have not yet seen clear evidence for 
the negative effective magnetic pressure instability in the present 
DNS. Once these various issues are better understood, it becomes 
timely to improve mean-field modelling. Obviously, all the mean- 
field models of the negative effective magnetic pressure instabil- 
ity have ignored realistic profiles of density and turbulent intensity. 
Also, the mean-field models should really be three-dimensional to 
include the possibility that magnetic structures break up along the 
direction of the mean field and form bipolar regions, as was seen in 
models of BKR. 
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